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Abstract. We construct a discrete form of Hamilton's Ricci flow (RF) equations for a 
d-dimensional piecewise flat simplicial geometry, S. These new algebraic equations are 
derived using the discrete formulation of Einstein's theory of general relativity known as 
Regge calculus. A Regge-Ricci flow (RRF) equation is naturally associated to each edge, 
£, of a simplicial lattice. In defining this equation, we find it convenient to utilize both 
the simplicial lattice S and its circumcentric dual lattice, 5*. In particular, the RRF 
equation associated to £ is naturally defined on a d-dimensional hybrid block connecting 
£ with its {d — l)-dimensional circumcentric dual cell, £* . We show that this equation is 
expressed as the proportionality between (1) the simplicial Ricci tensor, Rci, associated 
with the edge £ £ S, and (2) a certain volume weighted average of the fractional rate of 
change of the edges, \ £ £* , oi the circumcentric dual lattice, S* , that are in the dual 
of (.. The inherent orthogonality between elements of S and their duals in S* provide a 
simple geometric representation of Hamilton's RF equations. In this paper we utilize the 
well established theories of Regge calculus, or equivalently discrete exterior calculus, to 
construct these equations. We solve these equations for a few illustrative examples. 
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1. Ricci Flow on Simplicial Geometries 

Hamilton's Ricci flow (RF) technique has profoundly impacted the mathematical sciences 
and engineering fields [Tl[2l|3lll]. While RF has been used to solve the Poincare conjecture 
[5], the bulk of the applications of this technique have been limited to 2-dimensional 
geometries [6j . It is well established that the analysis of a geometry with complex topology 
is most naturally represented as a coordinate-free, piecewise-fiat simplicial lattice, e.g. 
as is routinely done using finite volume [T] , finite element [8] , Regge calculus [9l [10] and 
discrete exterior calculus [llj . In coordinate-based approaches, the difficulties in properly 
handling coordinate singularities, as well as the joining of multiple coordinate charts over 
the manifold, can be prohibitive. On the other hand, the utility of piecewise-flat simplicial 
geometries in analyzing the RF of 2-dimensional geometries is well established and proven 
to be effective |6l [12]. This utility should extend to three and higher dimensions, and 
preliminary work is underway in this direction |131 [T^ [TSj [TB J . 

Simplicial and combinatorial RF is well established in 2-dimensions; however, we are 
unaware of any formulation of the RF equations for a piecewise-flat simplicial geometry 
in dimensions greater than or equal to 3. It is the purpose of this manuscript to develop 
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such higher-dimensional Regge-Ricci flow (RRF) equations. These equations can be used 
in the analysis of topology and geometry, both (1) in the numerical analysis of RF and, 
perhaps more importantly, (2) as a mathematical tool to bound Ricci curvature in discrete 
geometries and to analyze and handle higher-dimensional RF singularities |17| [T8] . The 
topological taxonomy afforded by RF is richer when transitioning from 2 to 3-dimensions. 
In particular, the uniformization theorem says that any 2-geometry will evolve under RF 
to a constant curvature sphere, plane or hyperboloid, while in 3-dimensions the curvature 
and surface will diffuse into a connected sum of prime manifolds |19j . 

We consider in this paper the structure of Hamilton's Ricci flow equations for a piecewise- 
flat simplicial geometry, S, and its circumcentric dual lattice, S* . In Sec.|4]we use the sim- 
plicial geometry approach known as Regge calculus (RC) or its equivalent, discrete exterior 
calculus (DEC), to motivate our definition of a discrete representation of Hamilton's Ricci 
flow equation on each edge, A of the circumcentric dual lattice, S* . We will refer to as the 
RRF equations. Based on this definition, and using Thm. 2 and Cor. 2 on the double dual 
on piecewise-flat geometries |20j . we prove a corollary re-expressing these RRF equations 
as a corresponding equation for each edge, i, of the d-dimensional piecewise flat simplicial 
geometry, S. We refer to these equations as the simplicial RRF equations. 



Definition 1. We define the dual Regge-Ricci flow equation for any compact, piecewise- 
flat simplicial geometry, S, as an equation for each edge. A, in the circumcentric dual 
lattice, 5*, 



(1) 



/ Fractional rate of change \ 
of a circumcentric 
dual edge, 



\ 



X dt 



I 



I Ricci tensor associated \ 
to the circumcentric 
dual edge, 
Rex 



Corollary 1. The dual Regge-Ricci flow equation for any compact, piecewise-flat simplicial 
geometry, S, can be expressed as an equation for each edge, I ^ S, 



(2) 



/ The volume averaged \ 

fractional rate of change 

of the dual edge, 
^ldX\ 



.Xdt/t 



Ricci tensor associated to 
the simplicial edge, 
Rce 



Leading to a motivation for our definition and proof of the corollary, and by way of 
introducing the proper background and motivations for this work, we express Hamilton's 
RF equation in an orthogonal vielbein basis in Sec. [2] In Sec. [3] we provide an overview of 
the structure of the Riemann tensor (Rm), Ricci tensor (.Re) and scalar curvature (R) on a 
discrete piecewise-flat simplicial geometry. Following this in Sec. [4j we motivate Def. 1 and 
prove Cor. 1. In Sec. [5} we explicitly construct the 3-dimensional RRF equations and solve 
these analytically for two well known analytic examples, the 3-sphere and the 3-cylinder 
geometries. We conclude this manuscript in Sec. [6] by addressing a few important open 
problems. 
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2. Hamilton's Ricci Flow Equation 

The RF equation of a d-dimensional Riemannian geometry was first introduced by Hamil- 
ton 121], 

(3) %^ = -2^'^- 

where R^u are the covariant components of the Ricci tensor Rc, and g^u are the components 
of the metric tensor. It is our goal to reconstruct the RF equation of Hamilton on a 
piecewise-flat simplicial geometry, S, by using the formalism of RC [9j. To this end, and 
without loss of generality, consider a point p G {A4,g} where we can construct a time- 
independent orthogonal basis, {ea} , Va = 1, 2 . . . .d, i.e. 

(4) 

The contravariant components of the metric tensor in this basis, 

(5) 9 = eag'^'eh, 
will be, by construction, diagonal, 

(6) g'''' = diag{gn- 

In this specially-chosen basis (Eq. |4]) the time rate of change of the metric will be expressed 
solely in terms of the time derivatives of its d diagonal components. 

(7) g = eag''^ei, = eag'"'ea 

We will find it convenient for our purposes to express the RF equation (Eq. [3]) in its 
mixed-index form, i.e. as the components of a vector-valued 1-form equation in this basis 
at point p in the sense of Cartan geometry [22]. The components of these mixed-index RF 
equations take a simple form, 

(8) g''''gaa = -2R\. 

Since the metric at point p is diagonal then, g"""" = 1/gaa, and the RF equation takes a 
simple form, 

/ Fractional rate of \ _ (jaa _ _2^a 
V change of gaa ) ~ gaa ~ 



3. Curvature of a Piecewise-Flat ^-Dimensional Simplicial Geometry 

In this section we provide a brief review of the curvature on a piecewise-flat simplicial 
lattice with emphasis on its application to the RRF equation. Although a more complete 
and thorough description will appear in the literature [20], we felt that a review would 
make the derivation more complete and provide a slightly different emphasis needed for 
the derivation. In particular, we will outline the construction of the Rm, Rc and R on the 
simplicial lattice as well as its circumcentric dual lattice. We find it convenient to leverage 
the framework of RC [9l 1231 \2^ 1251 \26\ I27j . However, an essentially equivalent treatment of 
curvature on piecewise-flat simplicial manifolds is afforded by recent developments in DEC 
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[TT] . Both approaches, RC and DEC, appear to have independently encoded the geometry 
and curvature structure on discrete simphcial lattices in similar ways, and both utilize the 
circumcentric dual lattice in a fundamentally important way [20]. 

In the Euclidean form of RC the geometry of the interior of each of its d-dimensional 
simplicies is assumed to be flat Euclidean space, E'^. The geometry of each simplex in 
S is completely determined by the square of its ('^^^) = '^^'^2'"^^ edge lengths. Therefore, 
the entire lattice geometry and its curvature structure is completely determined by the 
squared length of all its edges, {ia}^ ^^'^ o^^y these edges. The set of all edge lengths of 
the simplicial geometry is the simplicial analogue of the metric, i.e. 

(10) g ^ {il}, 

in that any geometric quantity on the lattice, e.g. its lengths, areas, volumes, curvatures 
are defined solely in terms of these edge lengths and nothing else. 

RC was originally developed for general relativity as a piecewise-flat, coordinate-free 
simplicial representation of the 4-dimensional spacetime geometry. The geometry interior 
to each simplex is assumed to be a Minkowski flat spacetime. However, this theory is equally 
applicable for the curvature analysis of a d-dimensional simplicial spacelike geometry of 
interest here. Accordingly, we will assume that the geometry interior to each d-simplex 
is Euclidean flat space. In this transition, the formalism is otherwise unaltered from 4- 
dimensional spacetime geometries to purely spatial geometries. We will use this fact to 
implement RC and to construct the RRF equation in d-dimensions in this manuscript. 

The circumcentric dual lattice is fundamental to RC \23\ [2^ [29l |30j which is not the 
case for dual lattices of other kinds, e.g. the barycentric dual lattice. Consider a compact 
simplicial lattice S and its circumcentric dual S* . We note that for simplicity (and not out 
of necessity), we consider only compact d-dimensional simplicial geometries in this paper 
— proper boundary terms have been considered in RC |31J. Furthermore, we also assume 
here, for the same reasons, that 5 is a well-centered Delaunay lattice |28| and S* is its 
Voronoi lattice. The vertices, edges, triangles, . . . and d-simplexes in S will be denoted by 
ctq, o"!, (72, (Jd; respectively. The geometry of each simplex is determined uniquely by its 
(2) squared lengths; consequently the geometry of the entire simplicial lattice, {^^} G 5 
is a function of all the edge lengths. Similarly, we label the vertices, edges, polyhedra 
and polytopes of S* by ci^, cr^_i, (^d~2^ • • •' '^i^ ^0' respectively. The geometry of each of 
these dual cells owes its existence, and can be expressed as, functions of the squared edge 
lengths of S. The subscript on the simplicial elements is its dimension, e.g. fj^ G 5 is 
a A;-simplex; where the subscript on the circumcentric dual polytope is its co-dimension, 
e.g. cr^ G S* is the {d — /c) -dimensional circumcentric dual element of cjfc £ <S- Perhaps 
the most important property of the Voronoi and Delaunay lattices is that each element 
cjfc G 5 is perpendicular to its dual element cr^ G 5*. The d-volume formed by connecting 
the vertices of the fe-simplex, Ufc, to the vertices of the {d — A:) -dimensional polytope forms 
what we refer to as a d-dimensional hybrid polytope, Vkk*, of volume, 

, / , V f * fdY\ M *, k\{d-ky. , ,, 

(11) \Vkk*\ = {(JkWk) = / o-fcAo-fc= |cjfc||(Tfc| = 
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That this d-volume is a simple product of the volume of the dual-cell pairs reflects an 
important property of RC and is due to the orthogonality between the simplicial and 
circumcentric dual lattice. In addition, these d-hybrid volumes, when applied to all pairs 
of A;-dual cells {c^jO"/;}, provide a proper tiling of the discrete geometry. The curvature of 
the simplicial lattice has compact support on such hybrid volumes, e.g. the Rm is naturally 
defined on ¥22*, the Rc on Vn*, and Rc on Vbo* [HI [23 |26l |27] . We provide an illustration 
of the five hybrid volumes in the case of a 4-dimensional lattice in Fig. [T} 

In constructing the curvature tensors on a simplicial geometry, one can take various 
approaches. On a simplicial lattice, it is well known that the curvature is concentrated 
on the (d-2)-dimensional hinges h [9j. The form of this curvature is a conic singularity 
at the hinge. One approach used by researchers is to rigorously analyze this singularity 
structure, explicitly constructing a Dirac delta distribution for the curvature on the hinge 
h by embedding the hybrid block and its simplicial support into a {d+ l)-dimensional Eu- 
clidean space |23t|32j. This Dirac delta distribution is assumed to be uniformly distributed 
over the hinge, h, and thus integral curvature measures were derived. An alternative and 
equivalent description has recently been used [201 IMl ESI 123 123 EO]- First, we take the 
support for this curvature to be uniformly distributed over the entire hybrid block. Sec- 
ond, we use the Voronoi area, h* and only this area to define the support for the sectional 
curvature. Finally, we assume the geometry within each hybrid block, Vhh* > is an Einstein 
space [33]. These approaches yield curvatures that are proportional to each other, and can 
be transformed consistently to each other by taking suitable traces or identifying certain 
eigenvalues. 

In the remainder of this paper we assume the curvature within the y(d-2){d-2)* hybrid 
block can be treated as an Einstein space, and we define the Gaussian curvature by the 
parallel transport of vectors restricted to the perimeters of the dual Voronoi polygons, 
h* = cr^_2 where h* is the dual polygon to the ((i-2)-dimensional simplicial hinge, h. In 
this sense, the area of the dual Voronoi polygon, h* , becomes the unique weighting area 
associated to hinge h in order to define the Gaussian curvature. 



where €h is the deficit angle associated to hinge h, and h* is the area of the Voronoi 
plane dual to hinge, h. This V(^d-2){d-2)* hybrid cell is defined as the set of points in the 
lattice closer to the hinge h than to any other hinge and provides a reasonable support for 
this curvature. In the next three subsections we briefly reexamine the construction of the 
Riemann (Rm), Ricci (Rc) and scalar curvature (R) on a simplicial lattice beginning as a 
guide with their respective definitions from the continuum. 

For notational purposes, whenever we write a volume using a capitol V, it will be un- 
derstood that this represents a d-dimensional volume. It can represent a d-dimensional 
hybrid or reduced-hybrid volume. Its subscript will identify to which volume it corre- 
sponds. For example, the hybrid cell V(^d-2){d-2)* associated to hinge h will be written as 
Vh = V(d-2)(d-2)*- On the other hand, if we write, simply, h, in the correct context with no 
capitol, V, then it will mean the volume of h = \h\. Our notation follows similarly for all 




Gaussian curvature 
associated to hinge h 



= Kh 



Angle Vector Rotates 



Area Circumnavigated h* 
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Figure 1 . Each of the five rows in this table illustrates one of the five hy- 
brid 4-dimensional polytopes. The black edges (thick lines) in the diagram 
are in the simplicial Delaunay lattice, S while the red (thick grey) edges 
are in the circumcentric dual Voronoi lattice, S* . The blue (thin) edges are 
added in the right-most column and connect the vertices of the Delaunay 
lattice to the vertices of the Voronoi lattice thus forming a 4~dimensional 
polytope. Each of these 5 polytopes, when applied throughout all similar 
elements of S, provide a complete tessellation of the lattice spacetime. The 
symbols Vkk* in the left column completely characterize the hybrid polytope 
in the d-dimensional lattice. The first number, k, is the /s-dimensional sim- 
plicial element Ufc S 5, while the second starred number is its co-dimension k 
dual polytope in cr^ E 5*. In the center column, we show the A:-dimensional 
simplicial element in the black line, and its dual (d— A;) -dimensional element 
in S* is shown in red (grey). 
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other elements, e.g. I = \l\ for the length of an edge, h* = \h*\ for the area of the Voronoi 



polygon as in Eq. 12 



3.1. The Simplicial Rm. In this section we examine the assumption alluded to in the 
previous section that in RC the geometry's locally an Einstein space within each hybrid 
cell, and that the curvature is uniformly distributed within this block. To this end, one way 
we can construct the Rm on S is by examining the parallel transport of a vector around a 
closed loop. In particular, a vector, v := v"" Bb will ordinarily return rotated by an amount, 
eaSv"", when parallel transported around the boundary of the area element, A'^'^Cc /\edi 

(13) Sv'' = -2R\,dv'A'^^. 

This defines a constructive way to explore the components of the Rm through the rotation 
operator. 

(14) 11=^ e^ACb R^ e^Ae^ . 

rotation ^™ loop of 
bivector circum- 
navigation 

The curvature of a d-dimensional piecewise flat simplicial manifold is concentrated on 
its co-dimensional 2 simplicial hinges, h [9]. Given any hinge, h = ad-2 £ S, there will be 
at least three d-simplicies sharing hinge h. We assume here that there are n > 3 of these 
hinging on /i, and we label these simplicies {o"rf(z)} for i = 1,2, ... ,n. Let us consider a 
loop, C, within the interior of the collection of these n d-simplicies that encircles h only 
once. In addition, we consider a vector v at a point po along C. We place pq at the 
intersection of C and the ((i-l)-simplex , cr(i_i(l,2) that is shared by the two adjacent d- 
simplicies, (7(^(1) and C7d(2). The parallel transport of this vector within o"(i(l) is trivial since 
the geometry is flat. The angle v makes with the d{d + l)/2 edges of remains constant 
as it is parallel transported along C. Furthermore, as we continue to parallel transport 
the vector along C into (7(^(2) it is equally trivial. There is no curvature concentrated at 
their juncture ((Trf_i(l, 2)), i.e. the two adjacent d-simplicies can be embedded into the 
same flat Euclidean geometry. We can continue to parallel transport v along C through 
the remaining [n — 2) d-simplicies trivially; however, when we return to the starting point 
Pq on C we will find the vector rotated by the deficit angle e^, 

n 

(15) e/, :=27r-^^i, 

where the summation is over all a"rf(i) sharing hinge /i, and Qi is the hyperdihedral angle 
between the two {d — l)-dimensional simplicial faces a^-xii — l,i) and ad-i{i,i + 1) of 
simplex ad{i) sharing hinge h. Remarkably, this rotation bivector lies in the plane, h* = 02 
perpendicular to hinge h, and this rotation is independent of the of the orientation and area 
ofC. This captures the integrated curvature at hinge h. The curvature is a conic singularity 
the loop. The closed curve C can be shrunk to an arbitrarily small area, Aq, around 
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h without affecting the rotation in v under parallel transport. Therefore the sectional 
curvature can be made as large as one pleases by shrinking the area, 

(Conical Sectional \ / a i u. j. j. \ 
\ f Angle vector rotates \ en 

Curvature oj hmge = — = hm — = oo. 

^ g I \ Area circumnavigated J Ac->0 Aq 

Is there a natural area that we can assign to the hinge, /i? What area does RC use? 
To answer these questions, we identify three features of S that yield a unique curvature 
measure. First, the simplicial geometry is piecewise flat. Second, each element of simplicial 
geometry, Ufc G 5 is orthogonal to its circumcentric dual element, cr^ G 5*. Thirdly, it has 
been shown that the Regge equations as well as the Regge action are uniquely defined as 
distributed quantities over such circumcentric dual elements |24|, [25t [3l] . In his original 
paper on RC, Regge demonstrated that the curvature of a simplicial lattice is concentrated 
at each of the co-dimension 2 hinges, h = ad-2 G 5., and this curvature at each h \s a. 



conic singularity as illustrated by Eq. 16, However, in RC it appears that this curvature 
is defined to be distributed uniformly over the convex hull of the triangle hinge h and the 
polygon, /i* = G S* dual to h. Following the notation introduced in the end of last 
section, this area, h* G S* is central to the hybrid cell Vh ■= V22*, and this hybrid cell is 
the set of points in S closest to hinge h than to any other hinge. The area circumnavigated 
is, h* := \h*\ since RC defines the loop C to be the perimeter of the dual polygon h* . This 
area, h* and the deficit angle, e/j, are both functions of the of the squared edges, {£?} of S. 
This (1) dual area h* = \h*\, and (2) rotation, e/j, yield the hinge-based (or equivalently 
dual hinge-based) sectional curvature in RC, 

(17) K, = K,, = '^. 



In order to use the continuum equations, Eqns. 13 14 we need to construct an appropriate 
set of basis vectors within the Vh hybrid cell defined by the convex hull of h and h* . 
Following the development in Sec. [2| we will find it convenient to define a time independent 
d-dimensional frame of orthogonal bases 1-forms, {e°} for a = 1,2, . . . ,d, and corresponding 
basis vectors, {Ca} for a = 1, 2, . . . , d, that span a given hybrid polytope Vhh*, 



(18) ea = gabe^. 

We refer to the following basis as the ,^A-hybrid basis. We first identify an edge i € h 
and a dual edge \ € h* in Vh- if we connect the vertices of this A to the d vertices of the 
hinge, h this will form a reduced hybrid d-simplex that we will refer to as Vhx, i.e. the 
hybrid cell Vh restricted to A. The sum of all Vhxs for each X £ h* equals the volume of 
the non-simplicial hybrid cell. 



(19) Vh=J2^^^ 



hX- 

xeh* 



We define our orthogonal basis on each of these reduced hybrid cells, Vhx- In order to 
capture the orthogonality between S and S* within the hybrid cell Vh, we define basis 
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vectors of four kinds. 



(20) 



where, 



rrii 

t. 



rrihx, 



OiO, 



yi G {1,2,, 
A G h\ 



,(i-3}. 



01 is the circumcenter of edge £ £ h, 

02 is the circumcenter of any one of the triangles, 0"2, in h such that 02 3 i 

03 is the circumcenter of any one of the tetrahedra in h such that (T3 9 (T2 , 



• 0(i-3 is the circumcenter of one of two the (d-3)-dimensional simplexes in h such 
that Ud-z B (Td-i, 

• Od~2 is the circumcenter of the {d — 2)-dimensional simplicial hinge, h = ad-2, 
and, mix is the segment from d to the point on lambda, 0\, so that this line segment is 
perpendicular to A, i.e. it intersects A at the circumcenter of the ((i-l)-dimensional simplex, 
A*. The set of (d-2) mutually orthogonal vectors 



(21) 



span the hinge, h = (Td~2^ while remaining two vectors, 
(22) 

span the dual polyhedron, h* . The hinge can be expressed as a (d-2)-vector, 

1 



(23) 



h 



can be ej 
e£ Aei Ae2 A . . . Ae^^s) 



{d-2)\ 

where the ((i-2)-volume of the hinge is /i = \h\. In addition, the dual hinge can be equiv- 
alently expressed either as the dual of Eq. 23 or by using the last two basis vectors in 
Eq.[2| 

* ex A ehx 



(24) 



h* 



h*h . 



\ex A ehx\ 

The contravariant components of this flat-space metric, g, in this vector basis is, 

/ 



(25) 



9 



ab 



diag 



-2 -2 -2 -2 \-2 -2 



V 



hinge, h 



^In 2-dimensions there the hinge is a vertex and the space is a panned by only two orthogonal vectors, 
{\,fhh\\, i.e. there are no mi's. In three dimensions there are also no m^'s and the space is a panned by 
the three orthogonal vectors, {I, \,rhh\\- It is only in four and higher dimensions that we start generating 
the mj's. 
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A dual set of orthogonal basis 1-forms can be generated, 
(26) a;" = r7"''efe, 

and the covariant components of the metric are just, 



(27) Qab = diag 



( \ 



\ hinge, h h* j 

We are now in a position to construct the Rmy^ associated to hinge h by examining the 



rotation operator, T?., defined by Eq. 14 The sectional curvature, K}^, is naturally defined 
on the convex hull of h and h* , where the convex hull is just the hybrid volume, = V22* 
(as illustrated in Fig. [T] for the special case of 4-dimensions) . Furthermore, this is the only 
nonzero sectional curvature within the hybrid polytope V/j. It is the orthogonality between 
any element of the simplicial lattice, ct/j, and its corresponding circumcentric dual, cr^, that 
yields an Einstein space geometry within each hybrid cell, and that each of the d{d — 1) 
non-zero components of the Rmy^ is proportional to the sectional curvature -fC/i[25j. No 
matter what the orientation of the loop, S, we take encircling the hinge /i, the unit rotation 
bivector will always lie in the plane of h* and have magnitude e/i, 

(28) (.)^(S )=/»*£/,. 

Therefore, the only non- vanishing components of the simplicial Rra associated to hinge h 
are when the first two indices lie in the dual polygon, h* . While, it is true that in the basis 
chosen above the only non vanishing components of the Rm are in the Voronoi plane, /i*, 



(29) R^ 



h*' 



Nevertheless, this would not describe adequately the curvature within the hybrid cell. 
In particular, any basis chosen in general position would ordinarily (except for a set of 
measure zero) not have a single vector that lies in the hyperplane spanned by h. This 
would yield a fully-populated rotation matrix. In particular, this would be true for our 

specially constructed basis if we rotated each of the basis vectors |^ rhi,m2, ■ ■ ■ fhd — sj 

even slightly out of the plane of h. In this sense, the space within each hybrid cell, Vh , is 
a Einstein space [ 33J . Additionally, if we restrict the parallel transport of vectors around 
the perimeter of the Voronoi polygon, h* , then we have usual Einstein space expressions 
for the Riemann, Rruh, the Ricci, Rch and scalar R^ curvatures associated to hybrid cell 
Vh , 

(30) Rrrih = Kh = ^, 

(31) Rch = {d-l)Kh = {d-l)^, and 

(32) Rh = d{d-l)Kh = d{d-l)^. 
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Each of these Einstein space curvatures (Eqs. 30 32) associated with their hybrid ceh 
are uniformly distributed within this cell. This interpretation is consistent with other 
derivations [ 231 [32] with proper care in their normalization. Therefore, the Riemann tensor 
in Cartan language is, 

(33) Rruh = h* Rh h* = h* did - 1) Kh h* = h* did -l)^h*, 

n* 

where, h* = h* /h* is the unit 2-form generated on h* in hybrid volume, Vh- In the i\- 



hybrid basis this unit bivector is given by Eq. 24 

As an illustrative example, we can use these expressions for the Rnih and the hybrid 
proper volume, Vh, to derive the RC expression for the Hilbert action [SI |25]. In RC the 
geometric expressions are are often expressed as integrated quantities over the relevant d- 
dimensional hybrid or reduced-hybrid volume elements, Fig. [T] The Regge-Hilbert action, 
Ir-h, in general relativity, is one such integrated quantity. 

^^^^ ^ = ^1 ^^w. - ^E^'^^^ 



h 



(35) 



li^E (^^'^-'^m) \did-l) 



(36) = ^E^'^h. 

h 

3.2. The Simplicial Ricci Tensor (Re). In the last section we constructed the Riemann 
tensor and in this section we examine the simplicial Rc on S. This will involve the simplicial 
version of taking the trace of the simplicial Rm tensor in Eq. [30] We examine the Rc 
tensor here as it appears on the right hand side of the RF equation (Eq. [3]). This tensor 
was first introduced in ^23j, and later in context within the hybrid cells by [27J. We briefly 
reconstruct this tensor by utilizing the Einstein space geometry interior to Vh, and by 



using the £A-hybrid basis construction (Eq. 20) discussed in the previous section. This 



construction provides us with a geometric description of the trace operation on alattice 
geometry. 

In the continuum, the simplicial Rc is a 1-form valued 1-form, 
(37) Rc = e''Rabe'' = e''R\cbe\ 

This tensor is defined along some basis vector, e?, in the tangent space. This basis vector 
defines a common set of area bivectors for circumnavigation, i.e. S = (j'^^Cc A e^, one for 
each, c ^ b . In the previous section we defined the Rruh in RC on the dual Voronoi areas, 
h* = o"2_2 £ "S* of the hybrid cell, Vh- Each edge, A = cr^_| G 5*, is a common bounding 
edge of d{d — l)/2 of these dual polygons, h* . It is therefore natural to define a simplicial 
Ricci tensor, Rc\, along these dual edges, A. In this sense the trace of the Rm is a suitable 
sum over each of these Voronoi ares, h* sharing edge A. We illustrate this in Fig. [2] for the 
case of three dimensions. 
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Figure 2. We illustrate here the lattice geometry used to define the Rc 
tensor for a 3-dimensional simplicial lattice. The Rc tensor is naturally 
defined on an edge, A = cr* £ S* . Along this edge we define one of the 
unit triad vectors, ex of the ^A-basis, where ex = Xex as usual. In this 
illustration, we show the three dual polygons in red, G S* sharing edge A. 
Dual to each of these polygons is an edge of the triangle h = X* £ S. Along 
one of these edges, £ C h, we define the second of the triad vectors, eg = lei. 
Finally the normalized vector from the center of edge £ perpendicular to edge 
A, defines our final triad vector, eix = 'nitx^m- 



We construct the Ricci tensor, Rcx , associated to a dual edge A by properly tracing 
over each of the Ricci tensors, Rch, defined associated with each Voronoi polygon, /i*, 
containing dual edge A. The trace, or sum, will be over all h* restricted to A, or h*\x- 
The trace is accomplished using a weighted average as described for RC in [23 EHl and 
independently in discrete exterior calculus by [TT] . Here one considers Rcx to be uniformly 
distributed over the hybrid d-volume of the simplicial lattice associated with edge Vx and 
integrates this over this volume. This integrated tensor can be identified with the trace of 



the familiar hinge-based expression for simplicial Riemann tensor, Rnih in Eq. 31 with the 
inclusion of a proper weighting factor. In particular, 

(38) RcxVx= Y.^^'hVxh , 

hex* 

^ V ' 

simplicial form 
of the trace 
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where, 

(39) V, = rj AA* = ^ r 1 hhi, 

^ ^ hex* ^ ^ 

(40) yA/^ = (2) ^^A, and 

(41) Rc, = (d-l)g. 
Hence, 

(42) Rc, = ^^'^ ""'f' {Rch. 

In constructing Rcx we make use of reduced hybrid blocks, Vh*\, as well as reduced hinge 
volumes, h\ which are illustrated in the case of d = 4 in Fig. [3j In general, the reduced 
hybrid d-volume, Vh*x, is the fraction of the set of points interior to the hybrid d-volume, 
Vh* that are closest to the dual edge, A. This provides a decomposition of the hybrid cell 
into reduced hybrid cells. 



(43) Vh=Yl ^h. 



^hX, 

as illustrated in Fig. [2] Each reduced hybrid cell, Vh\, contain the circumcenters, Oj, that 
defined the £A-hybrid basis. Similarly, the reduced hinge is a ((i-2)-volume, A^, and is 
defined to be the set of points in the interior of the ((i-l)-dimensional simplex, cr^, that are 
closest to its bounding hinge, h, than to any other hinge in its boundary. 

The last two indices, c and b, of the Rm in the far right-hand side of Eq. 37 identify the 
components of a circumnavigated area, and in our case this will be the Voronoi polygons, 
h*^, sharing edge A. In the /lA-basis the two indices identifying this area element are b = X 



and c = mhx- As we saw the Rc tensor is associated with dual edge, A. Therefore, Eq. 37 



becomes, this fixed the index the index b = X and it is the tangent vector defining the 
second index of the Rc\. The other index, c, in Eq. 37 must be the component identifying 
the moment arm, ruhx, i.e. c = rrihx- There is a unique ruhx for each of the dual (2) 
Voronoi polygons, h*^, sharing edge A. We define Rcx by summing over these /lA-moment 
arms. Therefore, we would sum over all all polygons, h* , sharing edge. A, i.e. 

(44) Rc{...,\)=ea R'^cb {ex) = Ca R'^cx = ^ R'^'^n.x. 

Finally, since we showed that the simplicial Rm has compact support and uniformly dis- 
tributed over the dual hybrid polytopes, V(rf_2){d_2)* aiid the space is locally an Einstein 
space. The only non-zero vector component, a left in Eq. |44] is along A. Therefore, the 
simplicial Rcx is diagonal, 

(45) Rc{...,\) = Y, exR'^'^.x = ^ exR''\* = exY. 

h* h* h* 

"■\\ "-[A "|A 
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'D 



Figure 3. This figure illustrates the various geometric elements that were 
used to define the Rcx in 4-dimensions. The line segment A = UV is the 
edge in the circumcentric dual lattice, 5*, where we define the Ricci tensor. 
The tetrahedron A* = {ABCD} S 5 is perpendicular to A. If we connect 
the endpoints of A to the four vertices of A* the resulting 4-dimensional cell is 
the A-hybrid simplex of volume, Vx = {l/d)\\* . In 4-dimensions the hinges 
are triangles. We highlight one of the 4 triangles, h = {A^B^C} C 9A*, 
bounding A* by shading the triangle. The point O lies at the circumcenter 
of A*. It is the point equidistant from vertices, A, B, C and D. Point O also 
lies on the dual edge A. We refer to the tetrahedron formed by connecting 
O with triangle hinge, h = {A, B, C}, as the restricted 3-volume, Ajjj. There 
are four such reduced 3- volumes subdividing. A*. If we connect this shaded 
hinge, h, to edge A, the resulting simplex, will be the /lA-reduced hybrid 
block of volume, V\h = {l/d)Xh\. 



To complete the analysis of the right hand side of the RRF equation, we need also to 
define the Ricci tensor, on an edge i ^ S. In other words we need to express this 
£-based Ricci tensor associating it with an appropriate weighted sum of the A-based Ricci 



tensor equation (Eq. 42) onto the simplicial edge I. In essence, we are taking the double 
dual of Rc\ [20j . We will have to take the double dual of the left-hand side of the RRF 
equation in Sec. |4]in order to prove the Corollary 1. In particular, consider an edge ^ G 5 
as illustrated in Fig. |4j Similarly, one integrates the i?Q over the hybrid d-volume of the 
simplicial lattice associated with edge and equates this with the trace of the familiar 
hinge-based simplicial Ricci tensor {Rch) [25] . 

(46) RctV,= J]i?Cft% , 



simplicial form 
of the trace 
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where the sum is over all hinges, h sharing edge £, and 



(47) Vi = [f] ^r = ^:(;) heh*, 

(48) Vh^^ = ^£ = ^(2) ^"""^ 

(49) Rrrih = d(d-l)^. 

In order for us to obtain the right-hand side of Eq. [2]we make use of the following definition 
for the edge-based weighted average: 

(50) V /, • 

This yields an explicit expression for the Ricci tensor, Rci, associated to an edge, I, in the 
simplicial geometry, 

YjI^ Rch hi 

(51) Rce = ^ ^ — = {Rch)e, 



or equivalently, 

(52) Rci = (d - 1 



(eh) 



{h* 



3.3. The Simplicial R. The scalar curvature on the lattice is naturally defined at a 
vertex, v £ S [26j. The construction of the scalar curvature, i?^, on a simplicial lattice is 
similar to the derivation of Rc£ in the last section. The scalar curvature requires a proper 



weighted average over the hinged-based Riemann curvatures Rh of Eq. 32 As we observed 



in Sec. 3.1 , the Rrrih in RC is defined over the hinge-based hybrid volumes, the Ricci tensor 



we derived in Sec. |3.2| was naturally associated with an edge-based hybrid volume, Vi, and 
the scalar curvature R^ is naturally defined on the vertices, v G S, over a vertex-based 
reduced hybrid cell This trace requires a sum over a suitable reduced hybrid volume 
|26j, and in particular the volume, Vhv This reduced reduced hybrid volume is the set of 
points of the hinge-based hybrid cell, Vh that are closest to one of its vertices, v than to any 
other. We illustrate this reduced hybrid cell for the case of d = 4 in Fig. [5] In particular 
we find, 

(53) RvV,= ^RhVhv , 



simplicial form 
of the trace 
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Figure 4. In the left-most figure, we focus our attention on an edge (. = AB 
in a 4- dimensional simplicial lattice, 5, and illustrate the hinges h sharing 
this edge. In this particular case there are six such hinges, hi = ABHi, for 
i = 1, 2, . . . , 6. The orientation on edge i = AB induces an orientation on 
each of the six triangles sharing the edge. In particular, the triangle is a 
bi- vector in 5, where hi = \hi\h = \ABHi\ (^^pl^T^) • central figure 

we focus our attention on one of the six hinges, /ig, sharing edge i, and its 
Voronoi dual polygon /ig £ S*. This dual polygon, h^ is perpendicular to 
the hinge, /iq. The dotted line shown in the figure is simply the intersection 
of these orthogonal 2-dimensional surfaces. The orientation of the edge, 
i also induces an orientation on the bi- vector, h^ £ S* . In particular, 

h* = \h*\h* = h* (hiY = h* • The circumcenter of hinge h^ 

is labeled Oq . This point is equidistant from vertices A, B and Hq. The 
"/lA-moment arm" is the edge mh\ = OiOq = {1/2)1 cot (9) reaching from 
the center of edge i to the circumcenter of triangle h^. It is instrumental in 
the Cartan construction of the Regge-Einstein tensor |24j . The right -most 
figure shows the fraction of hinge h^ closer to edge i than to its other two 
edges. This edge-based triangle in S has area he = ABOq = 1/2 i'^ cot {9q), 
which is useful in defining the Rc in RC. 



where, 

(54) K = (^^~'v* = ^l^^~'Kh\ 

(55) Vh, = Yl (2) 

(56) Rh = d{d-l)^. 
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Figure 5. Here we show the reduced hybrid volume restricted to the 
vertex v = aQ £ S. This reduced hybrid d-dimensional polytope, h^, is 
formed by connecting the kite-shaped set of points in the triangle h G S 
closest to the vertex v to the dual area, h* £ S* . 

This yields an explicit expression for the Riemann scalar, Ry, associated to a vertex, v G S, 

Ylh, hv 

(57) Rv = = {Rh)v, 
or equivalently, 

(58) i?, = d(d-l)|^. 

3.4. Overview: the Rm, Rc and R. In the last three subsections, we constructed 
the Riemann, Ricci and scalar curvatures on a d-dimensional simplicial lattice, S. The 
coupling of elements of the Delaunay lattice to dual elements of the Voronoi lattice appear 
fundamental in their construction. The orthogonality between elements in cjfc G 5 and their 
circumcentric duals, G S* gave rise to relatively simple diagonalized expressions for the 
curvature expressions. This closed-form factorization of the homology and cohomology 
structures admits a unique interpretation of the curvatures in terms of hybrid blocks and 
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reduced hybrid blocks. These hybrid blocks are the product of Delaunay and dual Voronoi 
elements, 



(59) ^fc=(J ^kr^al 



We framed the curvatures by (1) uniformly distributing the curvature within their respec- 
tive hybrid cells, (2) restricting the domain of parallel transport to the scaffolding provided 
by the dual circumcentric lattice, and (3) identifying the interior of each hybrid block with 
an Einstein space geometry, where the mixed Ricci tensor is simply 



(60) Rc=^. 

a 



We showed in the last three sections that the Riemann, Ricci and scalar curvatures are 
naturally defined on the hinge-based, edge-based and vertex-based hybrid cells, respec- 
tively. In particular, the expressions for the curvatures that we will use, in part, to derive 
the right-hand side of the RRF equations in the next section, are as follows: 



(61) Rnih = {Rmh)h 



{(-h)h _ Qi, 

Wh ~ h* 



(62) Rci ={Rch)i ={d-l 



(63) i?„ ={Rh)v =did-l) 



and. 



Eh' he, (hl/h* 
(64) Rex = {Rc,)x = (d - 1) 't. , , , 



Additionally, we can project or assign any one of these curvatures to the vertices, edges 
and hinges of S as shown in Table 3.4 
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hinge, h 


edge, I 


vertex, v 


Rm 


{^h)h _ f-h 
{h*)h h* 




{<^h)t 




Rc 




^ ' {h*)e 




('^-1)1^ 


R 


did-l)l^,=d{d-l)^T^ 




d{d 1) S''^^ 



Table 1. The expressions for the simpUcial curvatures. The three boxed 
equations on the diagonal are particularly useful in the RRF equation con- 
struction. The expression for the Riemann scalar will be useful if we consider 
normalized Ricci flow. 



4. The Simplicial Ricci Flow Equation 



Definition 1. We define the dual Regge- Ricci flow equation for any compact, piecewise- 
flat simplicial geometry, S, as an equation for each edge, A, in the circumcentric dual 
lattice, S* , 

^ Fractional rate of change ^ / Ricci tensor associated \ 



(65) 



of a circumcentric 
dual edge, 



\ 



IdX 

X at 



/ 



V 



to the circumcentric 
dual edge, 
Rex 



/ 



Remark 1. We have shown here and in ^27J Ricci tensor is naturally associated with a dual 
edge, X £ S* . We consider this edge, A and the associated hybrid cell, V\\* . We also use the 



^A-hybrid orthogonal basis introduced in Sec. 3.1 Additionally, Sec. 3.2 and [27] showed 



that the mixed Ricci tensor is diagonal in A. Motivated by Hamilton's RF equation, Eq.|9j 
with mixed indices in an orthogonal basis, then the RRF equation associated with this 
dual edge is. 



(66) 



9XX 
gxx 



-2R\ 



-2 Rex. 



Using the gw component in the covariant metric for the £A-hybrid basis, Eq. 27 we obtain 
the dual RRF equation for A, 



(67) 



A 

A 



-Rex 



Since the geometry of 5 and S* is completely determined by the edges, {4} £ S, and 
nothing more, and since there are ordinarily many more dual edges. A, than there are 
simplicial edges, £, then the system of equations given by the dual RRF equation, Eq. 67 



will ordinarily be overdetermined and have no solution. In other words each A is a function 
of the Ps. Therefore we must map these dual RRF equation onto the simplicial edges. This 
can be done using an appropriate weighted average. 
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Corollary 1. The Hamilton Ricci flow equation for any compact, piecewise-flat simplicial 
manifold, S, can be expressed as an equation for each edge, I ^ S, 



(68) 



/ The volume averaged \ 
fractional rate of change 
of the dual edge, 



Ricci tensor associated to 
the simplicial edge, 
Rc£ 



Proof. We express the RRF equation by mapping the dual RRF equation (Eq. 67 ) onto an 
edge i, 



(69) 



(A/A), = {Rcx)i. 



It suffices then to consider the right-hand side of this equation. Using the double dual 
(Thm. 2 and Cor. 2 of |20] ) on Rc\, we note that 



(70) 



RciVe = ^RcxVex. 



We also note that the i — X reduced hybrid volume is symmetric in its indices, 

(71) Vex = Vxe, 

and we can expand the i hybrid volume into a sum of these reduced hybrid cells, 

(72) Ve = Y,Vn- 



Therefore, we see by the definition of the weighted average in Eq. 70 that the right-hand 
side of Eq. |69] is simply 



(73) 



{Rcx)e = Rce. 



Applying the edge-based weighted average to the LHS of the equation as we just did for 
the RHS of the equation yields the simplicial RRF equation, 



(74) 



IdX 

A at 



-Rcf 



□ 



5. The RRF equation in 3-Dimensions: Two Simple Analytic Applications 
In 3-dimensions the RRF equation obtains a relatively simple form since the edges, i, are 



the hinges, h. Consider £ ^ S and the associated RRF equation, Eq. 74 The right-hand 
side of this equation can expressed in terms of the deficit angle at £ and the Voronoi area, 

t, 



(75) 



RHS = Rce = 2 
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Figure 6. We illustrate here the restricted hybrid volume, Vu*^^, which is 
just the usual hybrid volume Vu* further restricted to one of the dual edges 
A. A is one of the dual edges in the boundary of dual polyhedron i* . 



The left-hand side is a weighted average of the fractional rate of change of the dual edges, 

A c r, 

(76) ^^^ = 4Ex^-' 

where Vi\ is the reduced hybrid volume shown in Fig. 6l and Vg = ^££* is the volume of the 
hybrid block. The mutual orthogonality of ttt-a, ^, and A yields a relatively simple formula 
for the volume of this reduced hybrid tetrahedron, 

(77) Vn = \iXmx. 
Therefore, the 3-dimensional RRF equation is simply 

(78) ^mAA = -4e,. 

Each dual edge, X £ £*, is a function of 9 simplicial edges of the two tetrahedra in S that 
share the triangle A*; consequently , 

(TO) ^.-Et^- 
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This yields a set of simplicial, nonlinear, coupled, first-order algebraic equations that 
describes the evolution of the edges of the simplicial geometry. In the next two subsections 



we will examine the solution of Eq. 78 for two simple models, the 3-sphere and the 3- 



cylinder, and demonstrate that they agree with their continuum analogues. 

5.1. 5 Cell, 16 Cell and 600 Cell Model of 3 D RRF. In this section we wih ex- 
amine the RF of a 3~sphere geometry by comparing the continuum evolution with the 
solutions of the RRF equations (Eq. [78]) . 

The metric describing the geometry of a 3-sphere in spherical coordinates is 

(80) ds"^ = Qij dx'dx^ = [dx^ + sin^{x)dn^) ; 

here, we let the radius of the sphere, a(t), depend on an external time parameter, t, and 
dVL^ = d'^9 + sin {9)d(fP' is the usual spherical line element. The only non-zero components 
of the mixed-component Ricci tensor for this geometry are 

(81) R^^ = R% = R^^ = {d-l)K = 2^. 

Therefore, Eq.( [9]) governing the RF for this geometry is a single equation for the radius, 

(82) = -4. 

Under RF this sphere collapses parabolically in time starting from an initial radius, = 
a(0) at time t = 0, 

(83) a'^it) =al- 4t. 

We now analyze the RRF equations for three simplicial representations of S"^. There are 
three regular polytopes in 4-dimensions ^35j whose boundaries are built of regular tetra- 
hedra and of which each is a homogeneous approximation to a 3-sphere. These polytopes 
are referred to as the simplex (5-cell model), the cross polytope (16-cell model) and the 
icosahedral pentagonal polytope (600-cell model), as depicted in Fig. [Tj The geometry 
of each of these three simplicial lattices are determined by a single length parameter, the 
length of the regular tetrahedron edge. We could let the edges of these lattice geometries 
vary independently; however, for the purpose of this paper we assume that they are all 
equal and, therefore, the polytopes remain regular during their evolution under RF. 



Here we explicitly construct and solve Eq. 78 for each of these three models and compare 



their evolution with the continuum. In these models the curvature is concentrated at the 
edges. The 5-cell has three regular tetrahedra sharing each edge, i, so that the dual Voronoi 
polyhedron, i* is a triangle. The 16-cell model has four regular tetrahedra sharing each 
edge, i, so that the dual Voronoi polyhedron, is a square. Finally, the 600-cell model 
has four regular tetrahedra sharing each edge, i, so that the dual Voronoi polyhedron, i* 
is a pentagon. The dihedral angle of a regular tetrahedron is given by 

(84) e = sec^^ (3) ^ 1.23096, 

so that the deficit angle ep at edge ip for each of these models is simply, 

(85) ep = 2TT-p9, 
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Figure 7. We illustrate here a Coxeter plane representation for each of the 
three regular 4-polytopes whose boundary are a simplicial approximation to 
the 3-sphere geometry. The left 5-cell polytope is a simplex whose boundary 
consists of 5 regular tetrahedra, three of each sharing each of its 10 edges. 
The middle figure is an illustration of the 16-cell cross polytope. It has 
16 regular tetrahedra in its boundary. Each of its 24 edges is the meeting 
place of 4 of these regular tetrahedrons. Finally the right-most figure shows 
a Coxeter plane representation of the 600-cell model. As its name suggests 
it has 600 regular tetrahedra in its 3-dimensional boundary. Each of its 720 
edges share 5 tetrahedra and in this sense this is often viewed as the higher- 
dimensional analogue of the icosahedron. This is the highest resolution 
homogeneous rendering of a 3~sphere geometry by a regular tetrahedral 
lattice. 



(87) ~ ^ ^P- 



where p = 3, 4 or 5 for the 5~cell, 16-cell and 600-cell model, respectively. In particular, 

5 - ceU p = 3 es f« 2.59031 
(86) 16 - cell p = 4 eie ~ 1.35935 

600 -ceU p = 5 6600 - 0.12839. 

The dual edges, Xp, for each these models are proportional to the edge length, 

6 

The moment arm, nip := mx^, in Eq. 78 is defined with respect to an edge, ip, and a dual 
edge, Xp, in the Voronoi polygon, i* . I'his line segment reaches from the center of the edge, 
£p, to the circumcenter of the triangle. A*. The moment arm, for each of our three lattices, 
is also proportional to the simplicial edge length, 

(88) mp = ^ep. 
The left-hand side of the RRF equation, 

(89) LHS=J2mpXp = ^^^{£l), 
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yields the RRF equation for each of these three models, 



How do we compare this with the continuum; in other words, how do we relate the radius, 
a, with the simplicial edge, ip? One way we can supply such a relation is by equating the 
surface 3~volume of the boundary of each polytope with the 3-volume of the 3-sphere of 
radius, a; i.e. 



(91) 



Volume of a 3-sphere \ _ f 3- Volume of the boundary of 
of radius a J y the p-cell with edge ip 



Here, Np = 5, 16 and 600 is the number of tetrahedra in the boundary of the 5-cell, 16-cell 
and 600-cell polytope; respectively. Therefore, we obtain the desired relation between the 
radius and the simplicial edge length, 



So that we may compare the RRF evolution with the continuum solution (Eq. 97), we 



substitute Eq. 92 into the RRF equation (Eq. 90), yielding 



We observe from this RRF equation that, as in the continuum Eq. 83 each of the three 
regular polytopes collapses parabolically in time starting from an initial effective radius, 
Oo = a(0) at time t = 0, 

2/3 



(94) a^{t) = 



6^/2 / V2Np 
p I 37r2 




In Table. 5.1 we show the percent deviation of the component of the Ricci tensor, R^^^^, 
from its continuum value Ry.y. = 2 as a function of the deficit angle for each of these 
three models. These results are consistent with the predicted second-order convergence 
of the simplicial Ricci tensor to its continuum value as a function of deficit angle. In 
particular, the logarithm of this fractional deviation for the 16 to the 600-cell model gives 
a convergence rate of approximately second order, {slope ~ 1.88). 

We are currently examining the solutions of the RRF equations on a 600-cell model 
where we allow all of its 720 edges to freely vary[36j. The RRF equations in this case 
form a sparsely-coupled set of 720 nonlinear first order differential equations for the 720 
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Model 


Deficit Angle 
( radians) 


Percent Deviation from the 
Continuum Value of R^^ = 2 


5-cell 


2.59031 


41.0% 


16-cell 


1.35935 


20.5% 


600-cell 


0.12839 


2.02% 



Table 2. This table demonstrates the convergence of the RRF equations 
to their continuum counterpart as a function of increasing mesh refinement. 
The 600-cell model reproduces the collapse rate in the square of the radius 
of the sphere within ~ 2%. The results from the 16-cell to the 600-cell 
model show a convergence rate in the deficit angle of order ~2 which is 
consistent with the non-commutativity of the finite rotations, (e^), as was 
argued in [24j. 



simplicial edge lengths. Our goal with this expanded model is threefold. First, we wish to 
reproduce the exponential uniformization proven by Hamilton for those perturbations of 
the edge lengths that keep the curvature positive. Second, we wish to analyze the behavior 
and stability of the RRF equations derived here. Thirdly, we can begin examining the 
role that modified RF and simplicial diffeomorphism freedom have on the stability and the 
accuracy of RRF. This model also has 1200 dual edges. A, so that we can examine more 
closely the behavior of the projection of the dual RRF equations onto the simplicial edges, 



5.2. Icosahedral 3 Cylinder Model of 3 D RRF. In this section we apply our RRF 
equations to a 3-dimensional cylinder of radius, R{t). The metric for a 3-cylinder is given 
by, 

(95) ds"^ = Qijdx'dx^ = dz^ + (t) d^^ , 

where z G {— cxd,cxd} is the axial coordinate, and dO?' = dO"^ + sin (0)^0^ is the usual 
spherical line element. We let the radius, R{t), of the 2-sphere cross sections of the 2- 
cylinder depend on an external time parameter, t. For this metric, the only non-zero 
components of the mixed-component Ricci tensor for this geometry are, 

(96) = = 

Therefore, Eq: [9] governing the RF for this geometry is a single equation for the radius, 

(97) = -2. 

Under RF this sphere collapses parabolically in time starting from an initial radius, Rq = 
a(0) at time t = 0, 



(98) R^{t) = Rl- 2t. 
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We now analyze the RRF equations, Eq. 78 for a particular lattice representation of the 
S-cylinder. In particular we represent each of the 2-dimensional spherical cross sections of 
the 3-cylinder by an icosahedron with edge length s{t). Adjacent icosahedra are connected 
by axial edges of length a{t) as illustrated in Fig. [sj Since adjacent icosahedra are joined 
together by 12 identical axial edges, a, then this 3-cylinder model consists of only one type 
of block - a regular triangular prism shown in Fig. [9] and Fig. 10 That this block is not a 
tetrahedron is not problematic. Ordinarily, the geometry of a non-simplicial block is not 
uniquely determine by its edge lengths, (i.e. it is non-rigid). However, the geometry of this 
prism is determined uniquely by its two edge lengths, s and a, because of its symmetry. 





a 




Figure 8. This illustrates the icosahedral 3-cylinder model we use to ap- 
proximate the RF of a 3-cylinder. In theis lattice, the cross section of the 
3-cylinder is approximated by a regular icosahedron with edge length, s. 
Adjacent icosahedra are joined together at their corresponding vertices by 
20 equal axial edges, a. In this model we examine the evolution under the 
RRF equations of these two edge lengths. 

The circumcentric dual lattice for this model is similar to a regular triangular prism 
lattice. The cross sections are regular dodecahedra of edge length 

(99) a = (\/3/3)s, 

and adjacent dodecahedra are connected to each other by axial edges 

(100) a = a, 

as shown in Fig. 10 , The dual lattice is composed of a single block - a regular pentagonal 
prism. The dodecahedrons are located halfway between each icosahedron, and therefore 
the adjacent dodecahedrons are connected by 20 equal axial edges of length a = a. Each 
a-edge is dual to a face of an icosahedron, while each a is dual to a rectangular face of the 
triangular prism. 

In this model, the only non-zero curvature is concentrated at the edges a. As we will 
see, there is no curvature along edge s . In order to calculate the curvature at edge a and 
edge s we need to calculate each of their deficit angles. We note that there are two dihedral 
angles of a regular triangular prism. The dihedral angle associated with the triangle edge, 
s, is a right angle, 

(101) = \, 
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while the dihedral angle associated to the axial edge, a, is equal to the internal angle of 
the regular triangle, 

(102) 0, = ^. 

Each axial edge a has five regular triangular prisms hinging on it, so the deficit angles for 
edges a are 

(103) e„ = 27r - 50a = ^, 

and since there are four prisms sharing each edge icosahedral edge s, the deficit angle 
vanishes, 

(104) es = 27T- AOs = 0. 
The only non-zero curvature is along the axial edges a. 



Figure 9. We show here a 3-dimensional illustration of the cross-sectional 
icosahedron and, on one of its 20 triangular faces, a regular triangular prism. 
This prism reaches from the icosahedron and will be attached to the corre- 
sponding triangular face of the adjacent icosahedron. Our lattice represen- 
tation of a 3-cylinder has only one polyhedral block, namely this regular 
triangular prism. 

The two edges, s and a, and their respective dual rectangle s* and dual pentagon a*, 
define three distinct moment arms. There are two moment arms associated with icosahedral 
edge, s. One of these connects the midpoint of s with the midpoint of the dual edge a, 

a 

(105) nisa = -, 
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Figure 10. We show on the left the triangular prism used to generate S. 
Due to the symmetry, its geometry is determined by the two edge lengths, 
s and o. On the right we show the polyhedron lattice block in the circum- 
centric dual lattice, S* . Its geometry is also fixed by assigning edge lengths 
to a and a. In the center we show the two polygons in S dual to the edges, 
a and a in 5*. 



and the other connects the midpoint of s with the midpoint of the dual edge a, 

(106) rusa = — = —s. 

2 b 

There are two of these moment arms on the rectangle s* . Finally, there are five identical 
moment arms associated with edge a. These moment arms reach from the center of edge 
a to the center of edge a; hence, 

(107) maa = ^. 

The left-hand side of the RRF equations, Eq. l78l for this model is given by equations. 



5\/3 

(108) LHSa = TTlA A = 5ma„ & = ss, 

(109) LHSs = rnx A = 2msa<y + 2msaa = — - {as + set) . 

Consequently, this yields the two RRF equations for this icosahedral 3-cylinder. 
(110) 

(111) 



dt 5^3 
a s 
a s 
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How do we compare this with the continuum; in other words, how do we relate the 
radius, r with the simphcial edge, s? One way we can supply such a relation is by equating 
the surface area of a 2-sphere of radius r with the surface area of an icosahedron with edge 
length s; i.e. 



(112) 



Area of a 2-sphere \ _ f Surface area of an 
of radius r J \ icosahedron with edge s 



or. 



Substituting this into Eq. |110 gives us the RRF equation in terms of an effective radius, 
r, tha 

(114) 



r, that is, serendipitously, the exact result as we found in the continuum, Eq. 97 



and change in the axial edge, Eq. [Ill will leave the cylinder geometry invariant. 

In 3-D and higher RF, singularities in the geometry may form. The type-1 and type- 
2 neck-pinch singularities examined in the dumbbell and javelin models by [SH\ [5^ can 
be examined in RRF by allowing the icosahedral edge lengths, s, and the axial edges, a, 
to vary spatially from one icosahedron to another. The evolution of the dumbbell and 
javelin model is currently being examined [37] using the RRF equations introduced in this 
manuscript. We have already observed neck pinching qualitatively in RRF, a more detailed 
comparison with the continuum models, analysis of conformal and diffeomorphic modes, 
and a stability analysis are currently underway [37j . 

6. Future Directions: Singularities, Surgery, Stability and the Simplicial 

DiFFEOMORPHISMS 

In this paper we have formulated a discrete form of the Hamilton RF equations for a 
piecewise flat simplicial geometry for any dimension, d. We imagine that this simplicial 
formulation of RF will provide new insights both numerically and, more importantly, in 
representing the singularity and soliton structure on piecewise flat simplicial lattices. It 
may be possible to use the RRF mathematics to explore singularity structure in three 
and higher dimensions in novel ways [30]. Perhaps there are minimal simplicial structures 
characterizing RF singularities of various types, in other words we can ask, "Given a RF 
singularity (Type-1, Type-2 etc.), what is the simplest simplicial lattice that exhibits this 
singularity?" 

Hamilton demonstrated that the RF equations will, under certain conditions, evolve to 
produce singularities (41] . He discussed some intuitive solutions for dumbbell geometries 
producing pinching singularities. For the RRF equations to produce a faithful representa- 
tion of Hamilton's RF, the lattice geometry will need to evolve so as to properly approach 
these pinching singularities and then to avoid them by surgery |4H [3] . The analysis of such 
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singularities using the RRF equations is one of the primary directions guiding our current 
research. We are in the process of analyzing the dumbbell model with 3-sphere topology to 
reproduce the recent results of [38l |39] using the RRF equations. Our preliminary results 
demonstrate neck pinching, although more analysis is necessary [37j . Once we demonstrate 
that the RRF equations reproduce quantitatively correct neck pinch singularity formation 
we will begin to examine a simplicial implementation of surgery. This will involve separat- 
ing the neck pinch into two separate lattices and then generate an appropriate simplicial 
cap for each of these, thereby continuing the RRF evolution to produce a connected sum of 
geometries. While we are only beginning in this direction, we believe it possible, at least in 
3-dimensions, to produce an automated surgery procedure for a lattice geometry in order 
to produce a complete RF evolution to a connected sum of Thurston-type geometries for 
a broad range of initial simplicial geometries and topologies [H] . 

We have begun our numerical implementation of the RRF equations in 3-dimensions 
on 3-sphere lattices [36]. In particular, we are currently analyzing the RRF of geometries 
with 3-sphere topology with everywhere positive Ricci curvature [36]. In this analysis ini- 
tial simplicial geometry is a perturbation of the lengths of each of the 720 edges of the 
600-cell regular polytope (Sec. 5.1 ) about their regular values so as to keep the Rc positive 
everywhere. Hamilton proved that non-uniformities in such geometries will converge to the 
regular 600-cell polytope exponentially fast [H] provided that the RRF faithfully repro- 
duces continuum RF. We are also examining Type-1 and Type -2 singularity formation in a 
simplicial representation of the dumbbell model j37| I42j . Our preliminary results for both 
models are qualitatively consistent uniformization; however, we need to carefully examine 
the detailed exponential fast uniformization predicted by Hamilton and to examine numer- 
ical behavior of the normalized RRF equations and their stability/instability. For small 
perturbation around the 600-cell model we examined the 720 eigenvalues of the Jacobian 
for the RRF equations and we find numerous positive real values indicating a stability issue 
|36j . Examining the stability of the RRF equations and providing suitable reformulations 
of these equations is the primary focus of our research. 

We are animated by three directions for this analysis. First, Perelman introduced a gauge 
fixed (conformal and diffeomorphism) modified RF evolution by introducing an entropy 
integral - an integral which produces uniformization via a gradient flow [5j. We ask, 
"What is the simplicial analogue of this for a piecewise flat simplicial lattice"? Second, an 
approximate simplicial diffeomorphism structure for a piecewise flat simplicial geometry has 
been understood [33] and successfully implemented in RC. This should enable us to examine 
a modified RRF equation using these approximate diffeomorphic degrees of freedom in 
order to improve the uniformization of the simplicial geometry (e.g. well-centered Delaunay 
lattices, maximally regular simplexes, etc.). Finally, there are many similarities between 
RF and issues involved in evolving black hole spacetimes that we believe can inform both 
fields of research. 

The numerical solution of Hamilton's RF equations in 3 and higher dimensions involve 
much the same issues as do the vacuum black hole solutions of the Einstein equations. 
In fact, the Einstein equations for a vacuum spacetime is just the Ricci tensor. First, 
Hamilton showed that one can encounter various singularities in RF which is also true for 
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the numerical evolution of black hole spacetimes. Secondly, our recent results suggest that 
the RRF equations have instabilities in certain cases [36l EZj as do the discretized version 
of the standard ADM equations in numerical relativity. Thirdly, Hamilton suggested a 
surgery procedure in order to integrate through the neck pinch singularities in much the 
same way as the moving puncture procedure or horizon excision procedure is used in 
numerical relativity to numerically evolve the spacetime geometry avoid the black hole 
singularity. Finally, in general relativity one utilizes the diffeomorphic degrees of freedom 
and conformal factors through the initial value data as well as a judicious choice of the 
lapse and shift functions during evolution so as to provide a stable solution of Einstein's 
equations that also avoids singularities (e.g. constant mean curvature slicing and minimal 
shear). In a similar vein, Perelman's entropy function utilizes diffeomorphic and conformal 
modes to provide a gradient flow for a modified RF equation. This suggests that experience 
gained in general relativity may inform similar issues in the numerical analysis of RF, just 
as experience with RF evolutions may offer insights into numerical relativity. We view this 
as an exciting avenue of research that may bridge two disciplines. 
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